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The information contained in the nucleon form factors of the energy momentum tensor (EMT) [1] allows to address, 
among others, questions like: How are the total momentum and angular momentum of the nucleon [2] shared among 
its constituents? Or, how are the strong forces experienced by its constituents distributed inside the nucleon [3]? 
. An answer is known only to the first of these questions. Deeply inelastic lepton nucleon scattering experiments 
show that at scales of several GeV 2 about half of the nucleon momentum (considered in the infinite-momentum frame) 
is carried by quarks, the other half by gluons. An appropriate tool allowing to deduce further information on the 
^ ' EMT form factors became available only during the last decade, namely generalized parton distribution functions [4] 
accessible in hard exclusive reactions [5-13], see [14-18] for reviews. 

The perspective to learn in this way about the EMT form factors motivates theoretical studies. The nucleon EMT 
form factors were studied in lattice QCD [19-22], chiral perturbation theory [23-25], and models such as the chiral 
. quark soliton model [26-28]. (Particular aspects of the EMT form factors were discussed in this model in [29-33].) In 
' this work we study the nucleon EMT form factors in the Skyrme model [34] . 
I/"") , The Skyrme model has important virtues. First, the model respects chiral symmetry which plays a major role in the 
■ description of the nucleon. Second, it provides a practical realization of the picture of baryons as solitons of effective 
meson fields motivated in QCD in the limit of a large number of colours N c [35]. Though A^ c = 3 in nature does not 
seem to be a large number, large- N c relations are often found to agree with phenomenology, although the approach 
has also its limitations [36]. Hereby it is also important to keep in mind that the chiral limit and the large- A^ c limit 
• do not commute [37]. Last not least, the simplicity of the Skyrme model makes it an appealing ground for numerous 
theoretical studies [38-49], see [50-52] for reviews. 

The purpose of the present study is threefold. First, the study of the EMT will shed some light on the model itself. 
I 1 Of particular interest in this context is the stability of the Skyrmion. Second, we calculate the nucleon EMT form 
Oh, factors in the Skyrme model which are of phenomcnological interest for the studies of the hard exclusive reactions. 
Third, we compare the Skyrme model results to the outcome from other chiral models, in particular, from the chiral 
quark soliton model [27, 28]. This comparison will help to recognize general features of the EMT form factors which 
are common in the large- N c soliton approach and therefore likely to be robust. 

For completeness we remark that the general chiral structure of the EMT was discussed in chiral perturbation 
theory and/or chiral models in [53-56]. Issues of pion EMT form factors in lattice QCD were addressed in [57, 58]. 

The note is organized as follows. In Sec. II we introduce the nucleon EMT form factors. In Sec. Ill we briefly review 
how the nucleon is described in the Skyrme model. In Sec. IV we derive the model expressions for the form factors. 
In Sec. V we discuss the densities of the static EMT, i.e. Fourier transforms of the form factors. In VI and VII we 
discuss the properties of the form factors and the numerical results, and conclude in Sec. VIII. The Appendices A-F 
contain remarks on the notation, technical details and proofs of the theoretical consistency of the approach. 
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II. FORM FACTORS OF THE ENERGY-MOMENTUM TENSOR 



The nucleon matrix element of the total symmetric EMT are characterized by three form factors as follows [1] 
(p>\T^(0)\p)=u(p>) 

where P = (p+p')/2, A = (p'-p) and t = A 2 . The normalizations (p'\p) = 2p°(27r) 3 <5( 3 )(p'-p) and u(p)u{p) = 2M N 
are used, and spin indices are suppressed for brevity. 

In QCD it is possible to define gauge invariant separate quark and gluon parts of the EMT having both decompo- 
sitions analog to (1) with form factors M®{t), M^(t), etc. In these decompositions, however, in addition the terms 
appear g M!/ c^ l,G (t) with cQ(t) — —c G (t) which account for non-conservation of the separate quark and gluon parts of 
the EMT, as only the total EMT is conserved. The quark (gluon) EMT form factors are related to Mellin moments 
of unpolarizcd quark (gluon) generalized parton distribution functions [2] which enter the description of certain hard 
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exclusive reactions [14-18]. In the following we shall focus on the total EMT. For our approach it is of importance 
that in the large- N c limit the nucleon EMT form factors behave as [16] 

M 2 (t) = 0(N°) , J(t)=0(N°), dx(f) - 0(N 2 ) . (2) 

These relations hold also separately for the quark- and gluon-part of the EMT. 

The form factors of the EMT in Eq. (1) can be interpreted [3] in analogy to the electromagnetic form factors [59] 
in the Breit frame characterized by A = 0. In this frame one can define the static EMT 

2V(r,s) = / 44 exp(iAr) (p 1 , S' |T^(0)|p, S) . (3) 

The initial and final polarization vectors of the nucleon S and S' arc defined such that in the respective rest-frame 
they are equal to (0, s) with the unit vector s denoting the quantization axis for the spin. The respective form factors 
are related to T^„(r, s) by (the following Eqs. (4, 5) hold also separately for quarks and gluons) 

J(t) + jJ'(t) = yd 3 re- irA £ ^ Sl r,T 0fe (r,s), (4) 
dr(t) + jd{(t) + ^d'!(t) = Jdhe-^T^v) (rV-^^) , (5) 

M 2 (t)-^(M 3 {t)-2J{t) + ld 1 (t)) = j±- Jd 3 re— A T 00 (r,s), (6) 

where the primes denote derivatives with respect to t. The form factors are renormalization scale independent and 
satisfy at t = the constraints 

M 2 (0) = -J- /d 3 rT o(r,s) = 1 , 

MN J 

J(0) = fd 3 re i ^s i r j T 0k (v, S ) = ^ 

di(0) = ~ J d 3 r T y (r) (rV - y = d 1 . (7) 

The first two constraints mean that in the rest frame the total energy of the nucleon is equal to its mass, and that its 
spin is 1/2. The value of d\ is not known a priori and must be determined experimentally. However, being a conserved 
quantity it is on equal footing with other characteristic nucleon properties like mass, magnetic moment, etc. [60, 61]. 

The components Tbo(r,s) and £ y rjTofc(r, s) describe the distribution of nucleon momentum and angular momen- 
tum, while the components Tjfc(r) characterize the spatial distribution of strong forces experienced by the partons 
inside the nucleon [3]. In fact, Tjj(r) is the static stress tensor which (for spin and 1/2 particles) is given by 

T ij (r) = a(r)^--^S i ^+p(r)5 ij . (8) 

Here p(r) describes the radial distribution of the pressure inside the hadron, while s(r) is related to the distribution 
of the shear forces [3]. Both are related due to the conservation of the EMT by the differential equation 

2 ds(r) 2s(r) dp(r) 

5 + + ir^ = o . (9) 

3 or r or 

Other important relations originating from conservation of EMT are the stability condition 

o 

dr r 2 p{r) = , (10) 



i) 



and integral relations for d\ in terms of p(r) and s(r) (see Appendix of Ref. [27] for further relations) 

di = - - M N [ d 3 r r 2 s(r) = - M N [ d 3 r r 2 p{r) . (11) 



3 J v ' 4 

It is worth to review briefly what is known about d\. For the pion chiral symmetry dictates df^ = —M^ [60, 61]. 
For the nucleon the large- N c limit predicts the flavour-dependence |d" + df \ ^> |d" — df \ [16] which is supported by 
lattice calculations [20-22] and results from the chiral quark soliton model [26-31]. Both lattice QCD and chiral quark 
soliton model yield a negative d® = d" + df. Also for nuclei df is found to be negative [3, 62]. In Ref. [27] it was 
conjectured on the basis of plausible physical arguments that the negative sign of di is dictated by stability criteria. 

The value of di is of interest for phenomenology, since it contributes to the beam charge asymmetry in deeply 
virtual Compton scattering [30] . HERMES data [9] indicate a df of negative sign though this observation depends to 
some extent on the employed model for the behaviour of generalized parton distributions in the small-x region [63] . 
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III. THE SKYRME MODEL 



In this Section we introduce briefly the Skyrme model. The Lagrangian of the Skyrme model is given by 

C= m M^t/cWt) + J-^ tr F [U\d»U), U\d„U)] [U\d»U), U\d v U)] + tv F (U - 2) (12) 

where U is the SU(2) pion field, is the pion mass, e is a dimcnsionless parameter, trF denotes the trace over the 
flavour-SU(2) matrices, and is the pion decay constant, with F^ = 186 MeV in nature. The parameters scale as 

F^O^J 2 ), e = 0(N~^ 2 ), m 7r =O(iV c ), (13) 

in the large N c limit, such that £ = 0(N C ) which guarantees the correct large- N c scaling of the theory. 

In the large- N c limit the chiral field U is static. Using the Skyrme (or "hedgehog") Ansatz, U = exp[iTe T F(r)] 
with r = |x| and e r = x/r, one obtains the soliton mass M so i = — Jd 3 x C as functional of the radial function F(r) 



M sol \F]=4ir drr 
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F 2 /2sin 2 F(r) 2 \ sin 2 F(r) /sin 2 F(r) , , 2 \ , m 2 F 2 



The soliton profile F(r) which minimizes (14) is determined by the following differential equation 



(14) 



r 2 2sin 2 F(r)\ ... s rF'Or) F'(r) 2 sin 2 F(r) sin2F(r) sin 2 F(r) sin 2F (r) mlr 2 sin F '(r) . . 

l F W + — + -ttt 2 - A ^r- 2 - A = (15) 



e 2 F 2 J w 2 e 2 F 2 4 e 2 F%r 2 

with the boundary conditions F(0) = tt and F(r) —> with r — > oo. These conditions ensure that the soliton field 
has unity winding number which is identified with the baryon number in the Skyrme model. 

In order to ascribe spin and isospin quantum numbers to the soliton one considers C/(x) — > A(t)U{~x)A~ 1 (t) in 
(12) with an arbitrary SU(2) matrix A = a$ + iar, introduces conjugate momenta 717, = dL/dai,, and quantizes the 
collective coordinates according to 7Tb — > —id /dab considering the constraint a 2 , + a 2 = 1. This yields the Hamiltonian 

H = M so} + ^ = M sol + ^ (16) 

where J 2 and I 2 are the squared spin and isospin operators, respectively. These operators act on the nucleon or A 
wave functions which can be expressed in terms of the collective coordinates cij, for example, the wave function for a 
proton with spin up is given by \p*) = (ai + ia^/ir, etc. [38]. Here & denotes the soliton moment of inertia 



2tt ( 

9 = — drr 2 sin 2 F(r) 




2 4F'(r) 2 4sin 2 F(r) 



F + 



(17) 



It is customary to treat F^ and e as free parameters. One way to fix them consists in adjusting F^ and e such that 
the masses of nucleon and A-resonance are reproduced [38, 39]. Hereby one first determines the function F(r) which 
minimizes M so \ and then projects on the proper quantum numbers. I.e. one considers the nucleon and A-resonance 
(and hyperons in the SU(3) version of the model) as different rotational excitations of the same object. 

In order to guarantee a consistent description of the EMT form factors we proceed as follows: We minimize M so \ 
and quantize collective coordinates as described above. In general the model expression for any observable A is 

A = A LO + A Iot , (18) 

where A^q appears in leading order of the large- N c expansion, while A Iot arises from soliton rotations and is suppressed 
by 1/N C with respect to A^q. For symmetry reasons A^q could vanish. If A^q ^ Owe shall neglect the contribution of 
A Iot (unless it happens to vanish anyway). In other words, we consider corrections due to soliton rotation if and only 
if the leading order gives a vanishing result. This step is consistent from a rigorous large- N c limit point of view, but it 
implies that we strictly speaking have to limit ourselves to the comparison of isoscalar and separately isovector flavour 
combinations of nucleon properties, see Table. I, which are in general of different order in the large- N c expansion. The 
results obtained below for the total EMT form factors refer, however, both to the proton and neutron. The reasons 
why other approaches are not favourable are discussed in App. C. 

The masses of the nucleon and A-resonance as they follow from (16) are of the type (18). (In this case -Arot is 
suppressed with respect to A^o by two orders in 1/N C .) In our approach Mm and Ma are degenerated and 

M N = M A = M so i = mm M so i[F] . (19) 
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Experimentally the nucleon and A masses differ by about 30% which indicates the typical accuracy we may expect 
(in the best case) for our results. The Skyrmc model is. in fact, observed to agree with phenomenology to within a 
similar accuracy [51]. From this point of view our approach is "within the accuracy of the model". 
In order to fix the parameters and e we choose to reproduce exactly the following observables 1 

M A + M N = 2M so i = 2171 McV , M A - M N = ^ = 293 McV . (20) 

This requires for m w = 138 MeV the following values 

F n = 131.3 MeV , e = 4.628 . (21) 

The experimental value of = 186 MeV is underestimated by 30% which is typical in this model. In order to gain 
intuition on the performance of the model with parameters fixed according to (20, 21) we consider several nucleon 
observables: isoscalar (1 = = proton + neutron) and isovector (1=1 = proton — neutron) electric and magnetic 
charge square radii (r^)/ and magnetic moments /z/, the axial coupling constant gA, and pion-nucleon sigma- 

tcrm (TirN- The latter is evaluated analogously to Ref. [39] yielding <j n n = 84 MeV, which agrees well with results 
from recent analyses [65] , see Table I, which indicate a larger value for a^N than earlier extractions [66] , and also with 
results from other models, see [67]. All results are summarized in Table I, and agree with phenomenology to within 
an accuracy of 30% with the exception of gA which, however, comes close to its estimated large- N c value [40]. 

It is of interest to consider nucleon observables in the chiral limit. For that it is necessary to specify what happens 
to the Skyrme model parameters F v and e as m n varies. In this work we shall keep these parameters fixed at their 
values in the physical situation in Eq. (21). Hereby we neglect the fact that in general F^ depends on m, [68, 69]. In 
this way we find that the nucleon (and A) mass are reduced in the chiral limit by 

M N\ p hysical point ~ M N\chiral limit =51 MeV, (22) 

which is within the range of the values considered in chiral perturbation theory. The isovector electric and magnetic 
charge square radii diverge in the chiral limit. Further results are shown in Table I. 

We remark that the description of the nucleon in the Skyrme model with the parameter fixing (20, 21) is — both 
for finite m n and in the chiral limit, see Table I — comparable to that in Refs. [38, 39]. 



TABLE I: The nucleon observables electric and magnetic charge square radii {r^i and magnetic moments and the axial 

coupling constant gA obtained in the Skyrme model using baryon masses according to Eqs. (20, 21), see text and Footnote 1. 



Quantity 


Skyrme model, results obtained here for 


Experiment 




m w = 


138 MeV 




-A/a + M N in MeV 


2068 


2171 (fixed) 


2171 


M A - A/jv in McV 


183 


293 (fixed) 


293 


F„ in McV 


131 


131 


186 




0.69 


0.60 


0.72 


/ 2 ,1/2 . t 


oo 


0.96 


0.88 


/ 2 \l/2 . , 


1.07 


0.85 


0.81 


/ 2 \l/2 . , 
( r m)l = l m fm 


oo 


0.96 


0.80 


M/=o in nuclear magnetons 


0.51 


0.66 


0.88 


/^/ = i in nuclear magnetons 


5.66 


3.70 


4.70 


M/=o/a i /=i 


0.09 


0.18 


0.19 


9A 


0.84 


0.73 


1.26 


a^N in McV 





84 


60-80 



Also other choices are possible. In this way, however, we preserve the "tradition" [38, 39] of using baryon masses for parameter fixing, 
and obtain a comparably satisfactory description of the nucleon, see Table I. We stress that in our approach Eq. (20) does not imply 
that Mjv and Ma are correctly reproduced. Instead, due to (19), one has Mjv = Ma = 1085 MeV which agrees to within 15% with 
the experimental values. Notice that the nucleon mass appears overestimated also in other soliton models. The origin of this problem is 
understood, and one way to solve it was discussed in [52] where it was demonstrated that quantum corrections due to meson fluctuations 
significantly reduce the soliton energy. A different way of understanding this problem, namely on the basis of notions and methods 
known from nuclear physics, was indicated in [64]. 
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IV. THE EMT FORM FACTORS IN THE SKYRME MODEL 



In the Skyrme model the canonical EMT, which — introducing the notation U 

dC 



is given by 



rpflU 



d v <t> a - C 



(23) 



is already symmetric. (Notice that in (23) one has to sum over a = 1, 2, 3; and </>o must not be considered as an 
independent field, since it is constrained by </>q + (j> 2 = 1.) For the respective components of the static EMT in Eq. (3) 
one obtains (c.f. App. C) 

F 2 /2sin 2 F(r) 



T 00 (r) 
T ofe (r,s) 



8 



+ F'(r) 2 ) + 



sin F(r) f sin 2 F(r) 



2e 2 r 2 

■2 



2F'{rf 



mlF 2 



1 — cos F(r) 



(s x r) 2 

while T^lr) is given by Eq. (8) with 
Fl (2 sin 2 F{r) 



pj(r) with pj(r) 



sin 2 F(r) 



12 6 



F'z 



AF'(r) 2 4sin 2 F(r) 



p{r) 
s(r) 



24 



,2 



+ F'(r) 2 + 



sin 2 F(r) ( sin 2 F(r 
6e 2 r 2 



2F'{r)'- 



ralFl 



1 — cos F(r 



F 2 sin 2 F(r) 
4 e 2 r 2 



F'(r) 2 - 



sin 2 F(r) 



(24) 
(25) 

(26) 
(27) 



In the large- N c limit Mjv ~ 0{N C ) while the components of the 4- momentum transfer behave as A ~ 0(N C - 1 ) and 
A' ~ 0(N®). Therefore \t\ <C Mfj. Considering this and the large- N c relations (2) one obtains from (4, 5, 6) 



Af 2 (t) 



t 



5M 2 N 



di(t) 



di(t) = 



1 

~Mn 
15M N 



d 3 r Too(r) j (»V=t) 



Jd 3 rp(r) 



J(t) = 3 / d 3 r pj(r) 



JiM-t) 



(28) 
(29) 
(30) 

with Too(r), p(r) and pj(r) as defined in Eqs. (24, 25, 26), and the Bessel functions jo(z) = and ji(z) = —j' (z). 

From the large- N c behaviour (13) of the parameters F n and e we see that Mn — 0(N C ) and = 0(N C ) such that 
M 2 (i) = 0(N°), J{t) = 0(N°) and di(<) = 0(N 2 ) holds in agreement with the general large-iV c result (2). 

Next we verify the internal consistency of the derived model expressions and demonstrate that they satisfy the 
general relations (7, 9, 10). For that we first take the limit t — > in Eq. (28) and use the fact that d\{t) is well-defined 
at t = 0, see below. By considering Eqs. (14, 19) we see that the first constraint in Eq. (7) is satisfied 

1 



M 2 (0) 



M 



N 



d J rT 00 (r) = 1 



(31) 



Similarly, taking t — > in Eq. (30) and comparing Eqs. (17, 25) we see that the second constraint in Eq. (7) is satisfied 

1 

2 ' 



J(0) - / d 3 r Pj {r) 



(32) 



To see that also the third constraint in Eq. (7) is satisfied it is convenient to prove first Eqs. (9, 10). To prove the 
stability condition (10) we rewrite r 2 p(r) as (for an alternative proof see App. B) 



r 2 p(r) 



d_ 

Or 



F'z 



rMr) + ^r^(r) 2 - ^ 



sin 4 F(r) F 2 ml 



G 



r 3 (l-cosF(r)) 



— ^- rF (r) x [ equations of motion ) 



(33) 



where "(equations of motion)" denotes the left-hand-side of Eq. (15) and vanishes. Thus, r 2 p(r) is a total derivative 
of a function which is zero at the boundaries of the integral in Eq. (10), i.e. J^°dr r 2 p(r) = is satisfied. Next, by 
inserting the expressions (26) and (27) into the differential equation (9) we see — in the notation of Eq. (33) — that 

F 2 

X _7T 

,.2 



2 ds(r) ( 2s(r) ( dp(r) 
3 



Or 



dr 



—^r F'(r) x (equations of motion ] = . 



(34) 



Now we are in the position to check the third constraint in Eq. (7) for di(t) at zero momentum transfer. Expanding 



the expression in (29) for small t, one obtains 



+ 0(t 2 ). The coefficient of the jr-term in this 



series is proportional to / °°dr r 2 p(r), i.e. identically zero. The next term in this series yields the expression for d\ in 
terms of p(r), cf. Eq. (11). The expression for d\ in terms of s(r) in (11) automatically follows from integrations by 
parts and the relation (9) which we have shown to be valid in the model, see Eq. (34). 
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V. PROPERTIES OF THE DENSITIES AND RESULTS 

It is instructive to consider first the Fourier transforms the nucleon EMT form factors (28-30), namely the "densities" 
Too(r), pj(f), p{r) and s(r) in Eqs. (24-27). Some of their properties can be derived analytically in the Skyrme model. 
For example, at small r the densities behave as 

T 00 (r)=T 00 (Q)-Ar 2 + O(r A ) (35) 

p(r) = p(0) -\Br 2 + C(r 4 ) (36) 

s(r) = \bt 2 + G(r 4 ) (37) 

5 

pj{r) = Cr 2 + 0(r 4 ) (38) 

which one concludes from solving the differential equation (15) iteratively at small r, see App. D. All constants Too(0), 
p(0) and A, B, C in (35-38) are defined to be positive numbers, see App. D for explicit expressions. 

Also the large-r behaviour of the densities can be determined analytically from that of the soliton profile function. 
Chiral symmetry and the variational equation (15) dictate the large distance behaviour of the latter to be as follows 

F(r) la '! c r ^ (1 + m„r) exp(-m T r) , (39) 

In the chiral limit the constant Rq, which is introduced here for later convenience, is related to the axial coupling 
constant qa according to 2 

R » = l§' («) 
Due to (39) the densities decay at large r exponentially for m T > 0, but exhibit power-like decays in the chiral limit 

T 00 (r) = 3i? 2 l + ... (41) 

,2 1 



p(r) = -R 2 -g + . . . (42) 
s(r) = 3i? 2 l + ... (43) 

Mr)= j§£ + ... (44) 



where the dots denote terms which are more strongly suppressed at large r. The results (41-44) allow to draw several 
interesting conclusions. Let us define the mean square radii of the energy and angular momentum density as 

(r |) = IggM , (45, 

_ jfcwi 

J d J r pj(r) 

We see that in the chiral limit: (r 2 E ) is finite, (r 2 ) diverges, di oc dr r 4 p(r) is finite, but the slope of di(t) at zero 
momentum transfer, ^(0) oc / °°dr r 6 p(r), diverges. 

The same small- and large-r behaviour of the densities was observed in the chiral quark soliton model in [27, 28]. 
(For the latter this is expected because (41-44) follow from the chiral behaviour of the soliton profile (39) equally 
respected in both models.) This indicates that the properties (35-38) and (41-44) are general features in the large- iV c 
chiral soliton approach and independent of dynamical details of underlying theories. 

Next we discuss the numerical results for the densities starting with the energy density T 00 (r) of static EMT (3). 
Fig. la shows TooM normalized with respect to the nucleon mass as function of r for the physical pion mass in 



2 The relation of to the asymptotics of the profile (39) is valid only in the chiral limit. For example, for the parameters used here 
Eq. (40) would imply for qa the value 0.76 to be compared with the correct result 0.73, see Table I, which follows from evaluating 
directly the axial current matrix element. We are grateful to the referee for drawing our attention to this fact. 
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comparison to the baryon density B(r) which coincides with the isoscalar electric charge distribution in the Skyrme 
model. Both curves in Fig. la have the same normalization and yield unity when integrated over the full space. The 
similarity of the curves for r < 1 fm means that nucleon mass and baryon charge are distributed similarly in this 
region. At large r, however, B(r) falls off much more rapidly than T 00 (r), in the chiral limit B(r) cx 1/r 9 vs. Eq. (41). 
Therefore the mean square radii, (r%) = (0.73 fm) 2 for the physical pion mass (and (0.94 fm) 2 in the chiral limit) in 
Table II, are much larger than the corresponding values of (r 2 ; ) r = o in Table I. The energy density in the center of 
the nucleon is 2.3GeVfm -3 , i.e. roughly 10 times the equilibrium density of nuclear matter. Figs, lb and c show the 
behaviour of Too(r) in the chiral limit compared to the physical case with emphasis on the small- and large-r regions. 

Fig. 2 shows the results for the angular momentum density p,j{r) which is related to the T ok components of the 
static EMT (3) according to Eq. (25). In the chiral limit p.j(r) is strongly decreased in the small-r region compared to 
the physical situation, see Fig. 2a. This is compensated by the slower decay of pj(r) in the large-r region, see Fig. 2b, 
which ensures in either case the correct normalization Jd 3 r p.j(r) = Sjy = 5, cf. Eq. (7). For the mean square radius 
(r 2 ) of the angular momentum density (46) we obtain (r 2 ) 1 / 2 = 0.96 fm at the physical point. In the chiral limit (r 2 ) 
diverges, see above. 

Finally we discuss the distributions of pressure and shear forces which characterize the T lk components of the static 
EMT according to (8). The pressure, see Fig. 3a, assumes at r = its global maximum with 3 p(0) = 0.48GeVfm -3 . 
It decreases with increasing r, initially according to Eq. (36), changes the sign at r$ = 0.64 fm till reaching its global 
minimum p m j n = -0.012GeVfm~ 3 at r min = 0.82 fm. For r > 

i"mm the pressure increases tending to zero but 
remaining always negative according to (42). The distribution of shear forces, shown in Fig. 3b, satisfies s(r) > 
and becomes zero only for r — > and r — > 00 according to Eqs. (37, 43). The shear forces exhibit one maximum at 
r max = 0.41 fm with s max = 0.15GeVfm -3 . The above-mentioned numbers refer to = 138 McV. In the chiral 
limit the distributions are qualitatively similar, see Fig. 3. 

The positive sign of the pressure for r < ro corresponds to repulsion, while the negative sign in the region r > tq 
means attraction. 4 In a mechanically stable object the repulsive forces in the inner region must precisely be balanced 
by the attractive forces in the outer region — in order to satisfy the stability criterion (10). In Sec. IV an analytical 
proof was given that the condition (10) is satisfied. It is instructive, however, to investigate in more detail how the 
stability is realized. For that we decompose the pressure in Eq. (26) into the parts which originate respectively from 
the 2- and 4-derivative term in the Lagrangian (12) — considering for simplicity the chiral limit where the mass term 
drops out. (The contribution of the 4-derivative term is labelled in Eq. (26) by the factor 1/e 2 .) 

The contribution of the 2-derivative term is negative for all r, see Fig. 4, which means that this term provides only 
attraction. This is an alternative way of stating that there is no soliton solution in the limit 1/e 2 — * in Eq. (12). 
The soliton collapses (shrinks to zero) in this limit as can be shown by means of scaling arguments, see e.g. [51]. 

The 4-derivative (Skyrme-) term — introduced, in fact, for the sake of stabilizing the soliton — is always positive, 
see Fig. 4, i.e. responsible for repulsion. For r < tq the repulsive forces from the 4-derivative term dominate over the 
attractive forces of the 2-derivative term. For r > r$ the situation is reversed. The two terms yield contributions to 
the stability condition (10) which are of opposite sign 

00 00 

J dr r 2 p(r)| 4 _ derivative term = - J dr r 2 p(r)| 2 _ derivative term = 13.7McV. (47) 


For m w ^ there is, in addition, a contribution of the mass term which is negative for all r and provides additional 
attraction. For the physical pion mass the decomposition analog to (47) reads 

{15.4MeV for i = 4-derivative term, 
-12.3 MeV for i = 2-derivative term, (48) 
— 3.1 McV for i = mass-term. 



VI. CHIRAL PROPERTIES OF THE FORM FACTORS 

Before we turn to the study of the chiral properties of the EMT form factors (28, 29, 30) it is important to remark 
that the limits N c — > 00 and — > do not commute [37]. The reason for that is the special role played by the 



This corresponds to 8 • 10 34 Newton/m 2 . One may gain some intuition on this number by considering that this is about an order of 
magnitude higher than the pressure in the center of a pulsar according to typical neutron star models. By multiplying p(0) with a 
typical hadronic area of, say, 1 fm 2 one obtains a force ~ O.SGeVfm" 1 which is about half of the QCD string tension. 
Studies in the chiral quark soliton model [27] suggest the following interpretation. The region r < ro is dominated by the "quark core" 
where Pauli principle is responsible for repulsion. The region r > ro is dominated by the pion cloud responsible for binding forces in 
the effective chiral theory, and thus attraction. Due to the lack of explicit quark degrees of freedom no such interpretation is possible 
in the Skyrme model but the results are in qualitative agreement with [27]. 
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FIG. 1: (a) The energy density in the nucleon normalized with respect to the nucleon mass, Tqq (r)/Afjv > as function of r in the Skyrmc 
model (solid line), and the baryon density B(r) for sake of comparison (dashed line). Both curves yield unity upon integration over the 
full space, (b) Too{r)/Mpf as function of r in the Skyrme model in the physical situation (solid line), and in the chiral limit (dashed line), 
(c) The same as (b) but for the normalized energy density 47rr 2 Too(f)/Mjv, i.e. here the curves yield unity upon integration over r. 
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FIG. 2: (a) The angular momentum density pj(r) normalized with respect to the nucleon spin Sjv = \ as function of r for the physical 
pion mass (solid line), and in the chiral limit (dashed line). The curves yield unity when integrated over the full space, (b) The same as 
(a) but for the normalized angular momentum density Airr 2 pj(r)/ Sn > i.e. here the curves yield unity upon integration over r. 
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FIG. 3: The distributions of (a) pressure p(r) and (b) the shear forces s(r) as functions of r in the Skyrmc model for the physical pion 
mass (solid line) and in the chiral limit (dashed line). 
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FIG. 4: The decomposition of the pressure in the chiral limit in Skyrme model. The graph shows the contributions to r 2 p(r) as functions 
of r due to the 2-dcrivative (dashed) and 4-derivativc term (dashed-dotted) in the Lagrangian (12), and the total result (solid line). 
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TABLE II: Different quantities related to the nucleon EMT form factors and their densities: the energy density in the center of the 
nucleon Tbo(O), the mean square radii (r^) and (rj) defined in Eqs. (45, 46), the pressure p(0) in the center of the nucleon, the position 
ro of the zero of the pressure defined as p(ro) = 0, the constant d\, the dipolc masses of the form factors M2(t), J(t) and d\(t) as defined 
in Eq. (55), and the mean square radius of the operator of the EMT trace {r? T ). In the chiral limit J(t) and di(t) have infinitely steep 
slopes at t = 0, see text, and dipolc fits do not provide useful approximations. For sake of comparison results from the chiral quark soliton 
model are shown [27]. 
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Skyrme model, results obtained here: 
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2.28 0.54 
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0.477 



0.70 
0.64 



-6.60 
-4.48 



1.00 

1.16 0.99 
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chiral quark soliton model, Ref. [27]: 




138 



1.54 0.79 
1.70 0.67 



1.32 



0.195 
0.232 



0.59 
0.57 



-3.46 
-2.35 



0.87 
0.91 



0.75 



0.65 



1.01 
0.81 



A-resonance. The mass difference Ma — Mjy = 0(N C T 1 ) vanishes for N c — > oo, and one has to consider the A on 
equal footing to the nucleon as an intermediate state in chiral loops. 

One consequence of the non-commutativity of these limits concerns the coefficients of the leading non-analytic 
(in the current quark masses) contributions (i.e. terms proportional to odd powers or logarithms of m,), which are 
model-independent — i.e. must follow from any consistent chiral approach. These coefficients depend on whether N c 
is kept finite or taken to infinity. For example the chiral expansion of the nucleon mass reads 



M N (m w ) = A + B ml - ki 



3gj 



with 



for N c 
for N r 



finite 
— > oo. 



(49) 



For scalar-isoscalar quantities the "discrepancy" is found to be always ki = 1 for finite N c vs. ki — 3 for large N c . 
For vector-isovector quantities a different factor ky appears which is ky = 1 for finite N c , and kv = f for N c — > oo. 
The factors fc, follow from soliton symmetries [71]. The form factors M2(t) and d\(t) are of scalar-isoscalar type. 

Scalar-isovector or vector-isoscalar quantities (J(t) is an example for the latter) appear at subleading order of the 
large- N c expansion and arise in the soliton approach from soliton rotations. For such quantities the non-commutativity 
of the limits N c — > oo and m T — > exhibits a more complicated pattern. A first example for that was encountered in 
the Skyrme model, where the isovector electric mean square radius diverges in the chiral limit as l/m ff [39] — and 
not as lnm ff [70]. The reason for that is that by limiting oneself to l/iV c -corrections due to soliton rotations only, 
one ignores pion-loop corrections which appear at the same order in 1/N C . The probably most advanced exploration 
of such corrections was given in [52] . Taking them here into account would go beyond the scope of the present study. 

Let us discuss first the chiral properties of the form factor d\(t). Its chiral expansion at t = reads, see App. E, 



d\ = di 



log 



2 M N 



16ttFJ 



■OK) 



(50) 



Here and in the following the superscript above some quantity denotes its value in the chiral limit. Thus, d\ is finite 
in the chiral limit. Noteworthy it receives a large leading non-analytic (linear in m 1T ) contribution. 

From the large-r behaviour of the distributions of pressure and shear forces in Eqs. (42, 43) we concluded that the 
derivative of d\(t) at zero momentum transfer diverges in the chiral limit. More precisely, we obtain, see App. E, 



<(o) = - 



0(m°) 



(51) 



For the chiral expansion of the mean square radius of the energy density, which is defined in (45) and which will 
be needed in the discussion of the behaviour of the form factor M 2 {t), we obtain, see App. E, 



(r 



81g 



16ttF?M n 



0{ml) 



(52) 
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Thus, (r E ) increases in the chiral limit — as all mean square radii do — but remains finite. The normalization of the 
form factor M2{t) at zero momentum transfer is fixed to unity, see Eq. (7). However, its derivative which is given in 
terms of (r 2 E ) and d\, see below Eq. (53), has a non-trivial chiral expansion, namely 

which follows from combining the results from Eqs. (50, 52). 

Also the normalization of J(t) at zero- momentum transfer is fixed, see Eq. (7), and what is of interest is e.g. J'(0) 
whose chiral behaviour follows from that of the mean square radius of the angular momentum density (46), namely 

J ' (0 » = i ir- "i> = ^ + °«'' ("> 

see App. E. We see that (r 2 j) diverges in the chiral limit which confirms our qualitative conclusion drawn from the 
large- r behaviour of the density pj(r) in Sec. V. In (54) one could eliminate the soliton moment of inertia in favour 
of the A-nucleon mass-splitting Ma — Mzv = 

The results for the chiral expansions in Eqs. (50-54) agree with the results obtained in the CQSM [27]. This is not 
surprizing. Both models describe the nucleon as chiral solitons whose structure is characterized by certain soliton 
profile functions. Those are determined by different dynamics, but what they have in common is the long-distance 
behaviour (39) uniquely fixed by chiral symmetry. (It is worthwhile remarking that it is possible to define a formal 
limit in which the Skyrmc model expressions follow from the chiral quark soliton model [72].) 

For d\{t) and M 2 (t) which appear in leading order in the large N c expansion our results agree with those obtained 
in chiral perturbation theory at finite N c [25] — up to the expected factor ki = 1 vs. 3, see Eq. (49). For J(i) which 
appears at subleading order of the large N c expansion we obtain J'(0) oc vs. lnm w at finite A c [25] which is 

also expected, see the discussion above. 

In the next Section we will discuss up to which values of the chiral expansions (50-54) provide useful approxi- 
mations for the exact model results. 



VII. RESULTS FOR THE FORM FACTORS 



Taking according to Eqs. (28-30) the Fourier transforms of the "densities" Tbo(r), p.j(r), p(r) discussed in Sec. V 
yields the nucleon EMT form factors. The large- N c approach can describe form factors for \t\ <€. M^. Here we restrict 
ourselves to < (— t) < lGcV 2 . In this region of t relativistic corrections were estimated not to exceed 30% [46] 
which is within the accuracy of the model. The results for the nucleon EMT form factors M 2 (t), J(t), di(t) are shown 
in Fig. 5. For the physical pion mass all EMT form factors can be well approximated to within an accuracy of (1-2)% 
by dipoles of the type 

F(Q) r 1.17 GeV for F(t) = M 2 (t) 

F ^ ~ 71 ,/\ ,2 x2 with M dip = { °- 99 GeV for F (t) = Jit) (55) 

U V^dipJ I 0.69 GeV for F(t) = di (t) 

for < (—t) < 0.8 GeV 2 . Fig. 5 shows the form factors also in the chiral limit. Here dipole approximations are not 
useful for J(t) and di(t) which exhibit infinitely steep slopes at zero momentum transfer. 

It is instructive to compare within the model the EMT form factors to the electromagnetic form factors — for 
example to the electric isocalar form factor G E =°(t) = G p E (t) + G n E {t) = G E (t) [45]. Fig. 6 shows G E (t), M 2 (t), J(t), 
di(t) normalized with respect to their values at t = where necessary. For the physical pion mass M 2 (t) and J{t) 
exhibit a similar i-dependence as G.e(t), while d\(t) exhibits a faster fall off, see Fig. 6a. The similarity of M 2 (t) and 
G_e(i) persists in the chiral limit, but d\{t) and J(t) have infinitely steep slopes at t = and fall off much faster with 
decreasing t than M 2 {t) and Gsif), see Fig. 6b. 

No principle fixes the normalization of the form factor d\{t) at zero- momentum transfer — in contrast to M2(t) 
and J(t). On the basis of stability requirements it was conjectured that d\ = di(0) must be negative [27], i.e. 

d x < . (56) 

In fact, the conservation of the EMT dictates J^dr r 2 p(r) = 0. For a mechanically stable object the pressure must 
be positive (repulsion) in the inner region r < r , and negative (attraction) in the outer region r > r , as we observe 
also in the Skyrme model, see Sec. V. This picture of the distribution of pressure forces, Fig. 7a, immediately dictates 
a negative sign for d\ oc J Q drr 4 p(r) = 0, Fig. 7b. Thus, the Skyrme model fully confirms the conclusion (56). 

Numerically we obtain d\ = —4.48 at the physical value of the pion mass for our choice of parameters, which 
confirms not only the sign but also the magnitude of results from the chiral quark soliton model [26-31]. 
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FIG. 5: The nucleoli EMT form factors M^ifjt J(t) and d\(t) as functions of t in the Skyrme model for the physical value of the pion 
mass (solid lines) and in the chiral limit (dashed lines). In the chiral limit J(t) and di(t) exhibit inhnitcly steep slopes at t = 0, see text. 




FIG. 6: Nuclcon form factors normalized with respect to their values at t = in the Skyrme model as functions of t: The EMT form 
factors M2(t) (solid line), J(t) (dashed line), d\(t) (dotted line), and the isoscalar electric form factor Gs(t) (dashed-dotted line) for the 
physical value of the pion mass (a), and in the chiral limit (b). 



r 2 p(r) in GeV fin" 1 r 4 p(r) x 57lM N in fin" 1 




FIG. 7: Relation of stability and the sign of di. (a) r 2 p(r) as function of r for the physical pion mass. Integrating the curve over r yields 
zero because the inncr-positivc-prcssurc-rcgion precisely cancells the outcr-ncgativc-prcssurc-rcgion. (b) r 4 p(r) multiplied by the constant 
jMj^Att as function of r. Integrating this curve over r yields di which apparently has a negative sign. 



It is also of interest to compare the chiral expansions in Eqs. (50-54) to exact model results which is done in Fig. 8. 
Figs. 8a, b show that the exact chiral behaviour of di and the mean square radius of the energy density in the model 
is dominated by the leading non-analytic contributions in Eqs. (50, 52) up to the physical value of the pion mass. 
Figs. 8c, d show that the slope of d±(t) at t = and the mean square radius of the angular momentum density 
diverge in the chiral limit as l/m v precisely as predicted by Eqs. (51, 54). (The smallest finite value of m n computed 
numerically in Figs. 8c, is about 1 MeV.) 

Finally, we turn our attention to the trace of the EMT. Due to the trace anomaly in QCD [73] the operator of the 
trace of the total EMT is given by 

ff = ^- F^F^ + (1 + 7m ) J2 ™ ^a . (57) 

" a 

Let F tT (t) denote the scalar form factor of the operator (57) which is, of course, not independent but a linear 
combination of the EMT form factors. The small-i expansion of F tI (t) = 1 +t(r^)/Q + 0(t 2 ) defines the mean square 
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FIG. 8: (a+b): The constant di and the mean square radius of the energy density as functions of m^. Solid lines: exact results. 

Dashed lines: Chiral expansions up to 0(m^) according to Eqs. (50, 52). (c+d): The numerical confirmations that ^(0) and (rj) diverge 
in the chiral limit as according to Eqs. (51, 54). 



radius (r| r ) which is related to (r 2 E ) and d\ as [27] 



J^rr^jr) 12 d, o 2 117g 2 



where the chiral expansion follows from Eqs. (50, 52). Numerically we obtain 



; 2 \ _ | 1.47 fm 2 in the chiral limit, /rq\ 
' 0.93 fm 2 for m w = 138 MeV. 



In the chiral limit (r 2 r ) is the mean square radius of the gluonic operator F^F^. The large value obtained in (59) 
is rather interesting, in particular, if we confront it with the mean square radius of the traceless part of the EMT 
estimated by means of QCD sum rules to be about 0.1 fm 2 [74]. A possible explanation why the radii of the trace 
and the traceless part of the EMT differ so much could be provided in the instanton vacuum model [75] , see [27] for 
details. 



VIII. CONCLUSIONS 



In this work we presented a study of the form factors of the energy momentum tensor of the nucleon in the Skyrme 
model. We have shown the theoretical consistency of the approach and provided explicit proofs that the EMT form 
factors in the Skyrme model satisfy all general requirements, provided one evaluates the respective expressions with 
that profile function which minimizes the actual expression for the nucleon mass. 

We derived the chiral expansions of the form factors in the Skyrme model and found them in agreement with 
results from chiral perturbation theory [24, 25] considering the non-commutativity of the limits m T — > and N c — > oo 
[37, 71]. The Skyrme model also confirms the strong m^-dependence of the constant d\ observed in the chiral quark 
soliton model [27, 28]. This is of practical interest in the context of extrapolating lattice QCD data to the physical 
value of the pion mass, see [28]. 

The numerical results for the form factors are in good qualitative agreement with results from the chiral quark 
soliton model [27]. Also the Skyrme model yields a negative constant d\ and confirms the conclusions drawn in 
[27] that di must be negative for a mechanically stable object. Numerically we find for the constant d\ = —4.88 in 
agreement with estimates from the chiral quark soliton model [26-31]. 

Results from both models suggest that all form factors can be well approximated by dipole fits for < (— t) < 
1 GeV 2 . M%{t) and J(t) exhibit a similar i-dependence as the isoscalar electric form factor, while d\{t) shows a much 
faster fall-off with decreasing t. Both the Skyrme and chiral quark soliton model [27] predict for d\(i) a dipole mass 
of 0.7 GeV. The predictions for the ^-dependence of d\(t) could contribute to an explanation of the t-dependence of 
the beam-charge asymmetry observed at HERMES [12]. Of particular interest is also the predicted i-dependence of 
the form factor J(t) which is of relevance for extracting from data the quark angular momentum contribution to the 
nucleon spin [76]. From the experience with the description of other observables in the Skyrme model we expect our 
predictions to hold to within an accuracy of (30-40)% though, of course, as always data will have the final say. 

In this work we used the simplest version of the Skyrme model [34]. It would be interesting to study the nucleon 
EMT form factors in extended versions of the Skyrme model which include e.g. vector- mesons [47-49], in order to see 
how these degrees of freedom contribute to the stability of the soliton and to the constant d\ . 

Acknowledgements We thank Pavel Pobylitsa and Maxim Polyakov for fruitful discussions and valuable comments. 
This research is part of the EU integrated infrastructure initiative hadron physics project under contract number RII3- 
CT-2004-506078, the Transregio/SFB Bonn-Bochum-Giessen, and the COSY-Julich project, and partially supported 
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APPENDIX A: ALTERNATIVE DEFINITION OF FORM FACTORS 

By means of the Gordon identity 2Mjs[u'j a u = u'(ia aK A K + 2P a )u one can rewrite (1) as (see e.g. in Ref. [2]) 



(j/\f$ G (0)\p) = u(p') 



A®- G {t) 



AM 



N 



C QG {t) A ^ A \ r 9 ^ A2 ±c(t) 9tw 



M 



N 



(Al) 



where A a (t) = M£(t), A a (t) + B a (t) = 2J a {t), C a (t) = with a = Q, G. In this notation the constraints (7) 

read ^4^(0) + ^4 G (0) = 1 and B®(0) + B G (Q) = meaning that the total nucleon gravitomagnetic moment vanishes. 

APPENDIX B: ALTERNATIVE PROOF OF THE STABILITY CONDITION 

The stability condition (10) can be proven alternatively in the following way. Let F(r) denote the profile function 
which minimizes the soliton energy (14), i.e. which satisfies (15), and substitute r — ► r' = Xr in (14). This yields 



M sol [A] = \- 1 E 2 + XE 4 + \~ E, 



(Bl) 



where E2, E4 and E m are the contributions of respectively the 2-derivative, 4-dcrivative and mass term in (14). 
M so i[A] in (Bl) has a minimum at A = 1, i.e. we have the conditions 



dM sol [A] 



<9A 

5 2 M sol [A] 



d\ 2 



A=l 



— E2 + E4 — SE m — , 
2E 2 + 12E rn > , 



(B2) 
(B3) 



which holds since the E\ (i = 2, 4, m) are positive. The above considerations are well known, the relation (B2) is 
sometimes referred to as the "virial theorem" in the Skyrme model [50, 51]. For our purposes it is sufficient to observe 
that the virial theorem is nothing but the stability condition 



j<lrr*p{r) = ±- 



1 9Af sol [A] 



dX 



= . 



(B4) 



A=l 



APPENDIX C: ROTATIONAL CORRECTIONS FOR T 00 (r), p(r) AND s(r) 



The expression for a quantity in the Skyrme model is in general of the type (18). In this work we have chosen to 
consider rotational corrections to a quantity if and only if it receives no contribution from the leading order of the 
large- N c expansion. This Appendix aims at making the arguments in favour of this procedure more clear. 

In our approach we have, in particular, neglected rotational corrections to the densities Too(r), p(r) and s(r). What 
would have happened if one kept these rotational corrections? Let us denote the densities with included rotational 
corrections by Tbo(r), p( r ) an d s(r). The corresponding model expressions read 



Too(r) 
p\r) 
s(r) 



T 00 (r) + 



S(S + 1 ) 
129 2 



sin 2 F(r) 



P(r) + sin 2 F(r) 



F 



F: 



4F'(r) 2 


Asm 2 F(r)~ 


e 2 


2 2 
err 


AF'(r) 2 


4 sin 2 F(r)~ 


3e 2 


3e 2 r 2 


SF'(r) 2 


4sin 2 F(r)" 


e 2 


e 2 r 2 



5(5 + 1) 
T 00 (r) + — pj{r) 



(CI) 
(C2) 
(C3) 



with Tqq(t), p{r) and s(r) as given by Eqs. (24, 26, 27), i.e. the corresponding densities without rotational corrections 
Note that Jd 3 x f 00 (r) 
considered in literature 



Note that Jd 3 xfoo(f) = M so \ + with M so i defined in (14). At this point several procedures have been 
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A. One procedure, used e.g. in [38, 39], consists in minimizing 
only the soliton energy M so \ in Eq. (14), but keeping the respective 
rotational corrections understood as small perturbations. Then 
e.g. the nuclcon or A mass are obtained by adding up the mini- 
mized soliton energy M so i and the rotational energy 3 20 ■ This 
procedure strictly speaking does not yield a stable soliton. This 
can be demonstrated e.g. by evaluating the pressure p(r) with the 
profile function which minimizes M so \. For the parameters used 
in [38] we obtain the result shown in Fig. 9. Clearly, the p(r) 
evaluated in this way is always positive and violates the stability 
condition (10). This is not surprizing because from the proof of the 
stability condition in the Skyrmc model given in App. B it is clear 
that (10) will be satisfied if and only if one evaluates the pressure 
with that profile function which minimizes the full energy. 

B. One way to resolve the problems encountered above consists 
in minimizing M so \ = Jd 3 x Too( r )- The profile function which 
does that has to satisfy the integro-differential equation (since O = 
G[F], see Eq. (17)) 



r 2 p(r) in GeV fm" 1 




-1 - 
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1 .... 1 . . 
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FIG. 9: r 2 p(r) in the Skyrme model as function of r. 
Dashed line: Resuit obtained following the procedure 
(A), see text, used e.g. in Ref. [38]. Solid line: Result 
obtained following the procedure (C) used in this work. 
Both curves refer to the chiral limit. 



left-hand side of Eq. (15) 



5(5 + 1) 
68 2 

8rsin 2 F(r)F'(r) r 2 
1?FZ + 1> 



4r 2 sin 2 F(r) F"(r) 2r 2 sin2.F(r) F'{rf 



sin2F(r) 



e 2 F 2 

4sin 2 F(r) sin2F(r) 



(C4) 



= 



with the same boundary conditions as for Eq. (15). Now one obtains a consistent description of the densities. E.g. 
the distributions of pressure and shear forces with rotational corrections satisfy the relation (9), because 



2 da(r) 

3 dr 



2s{r) 



or r^ 



equations of motion with rotational corrections J 



0. 



(C5) 



where "equations of motion with rotational corrections" denotes the left-hand-side of Eq. (C4). Also the stability 
condition (10) holds with included rotational corrections which can be seen by rewriting 



sm'F(r) , ^>i r3(1 _ cos ^ (r)) 



r z p(r) = |-(r 3 p(r) + §r a F'(r 



dr 



12 



5(5+1) 
69 2 



3 



3e 2 r 
4r 3 
3^ 



^-r 3 sin 2 F(r) + ^ sin 2 F(r) F'(r^ 2 



— ^-rF'(r) x (equations of motion with rotational corrections J . 



(C6) 



Alternatively, one may use the scaling trick discussed in App. B. For that we notice that in addition to the Ei in 
Eq. (Bl) there are two contributions due to rotational corrections, namely 



E T . 



S(S + l) r2 



F$ / d 3 x sin 2 F{r) , E mtj2 



s(s + i) r Ai 



J d 3 x sin 2 F(r) \F' ' {rf + 



sin 2 F(r) 



129 2 T ./ w ' ' 3e 2 2 

such that 

Msoi[A] = A _1 i?2 + Ai? 4 + A - E m + X~ 3 E Iot ^ + A _1 i? roti2 

Also in this case the conditions analog to (B2) hold and, in particular, we have 

00 

1 5M sol [A] 



dr r p(r) 



12tt dX 



= 



(C7) 



(C8) 



(C9) 



A=l 



Thus, also with included rotational corrections one obtains a consistent description of the form factors of the EMT. 
The problem is, however, that (C4) has no solution in the chiral limit, while the solutions which exist for sufficiently 
large yield results in conflict with chiral symmetry [43, 44], see also the review [51]. 

The minimization of M so \ = M so \ H — 2e would be unsatisfactory also from the point of view of large- N c counting, 
as no contribution of 0{N®) is included. Such terms could in principle be taken into account by adopting the 
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Peierls-Yoccoz projection technique in non-relativistic many body theory [77-79]. If one assumes the "small overlap 
approximation" the energy is similar to the one obtained approximately by pion corrections to the Skyrmc model 
[80]. In the chiral quark soliton model, which is in this respect similar to the Skyrme model, this has been discussed 
in Ref. [64]. 

C. A consistent description of the nucleoli EMT which respects chiral symmetry is obtained by adopting the 
procedure used in this work. After quantizing the soliton one neglects rotational corrections with respect to the 
leading order contribution unless the latter is zero. I.e. rotational corrections are considered if and only if the leading 
order contribution happens to vanish for symmetry reasons. 

In this way one strictly speaking cannot describe the proton and neutron separately, but has to limit oneself to 
isoscalar and isovector quantities which are in general of different order in the large- N c expansion. This is what we 
did in this work, see e.g. Table. I. However, the results obtained here for the total EMT form factors refer both to 
the proton and neutron if one neglects isospin violation (and electromagnetic corrections). 



APPENDIX D: BEHAVIOUR AT SMALL r 



Solving the Euler-Lagrange equation (15) iteratively at small r one obtains for the profile function the expansion 

Ok 

fc! 



F(r)=Y,lj( eF ^ k ( D1 ) 



fc 

where = it and = for even k > 2. With the definition = m- K /{eF 7r ) the afc for odd k are given by 

cti = at 

4a 3 (l + 2a 2 ) + 3am 2 
° 3 ~ 5(1 + 8a 2 ) 

24a 5 (5 + 32a 2 + 88a 5 + 448a 6 ) + 3a 3 m 2 (35 - 16a 2 + 832a 4 ) + 3am 4 . (5 - 56a 2 ) . . 

° 5 = 35(1 + 8a 2 ) 2 ' (D2) 

etc. Notice that a = a(m ff ). Numerically we find a = —1.00376 for fh n = 0, and a = —1.10752 for m, = 0.22715 
(which corresponds to = 138 MeV for our choice of parameters). In this way (15) can be solved in principle as a 
differential equation to any finite order in r. However, this iterative procedure can only be used after (15) is solved as 
a boundary value problem (with the constraints F (0) = ir and F(r) — > as r — > oo) and a = F' (0)/(eF„) is known. 

Inserting (Dl, D2) into the model expressions (24-27) for the densities Too(r), pj(r), p(r) and s(r) one obtains the 
results in (35-38) with the positive constants 

Too(0) = ^(3a 2 + 12a 4 + 4m 2 ) , A = — ^ a 2 (a 2 + 4a 4 + m 2 ) (D3) 

,(0)^(4a 4 -a 2 -4„ 2 ), B= Z* Z * ± *3 ( D4 ) 

_ e 2 F 4 a 2 (8a 2 + l) 
6 _ ~4 36 ■ (D5) 

From (D4) it is not obvious that p(0) and B are positive. In order to check that this is really the case, we note first 
that if p(0) > 0, then also B oc (8a 4 — a 2 — 2m 2 ) > (8a 4 — 2a 2 — 8m 2 ) oc p(0) > 0. Thus it is sufficient to demonstrate 
that p(0) is positive which can be rewritten as 

Tfl 2 1 

p(0) >0 o a 2 > ^ + - . (D6) 
a z 4 

By numerically solving the boundary value problem for different m, < 0.5, which for our choice of parameters in (21) 
corresponds to < 300 MeV, we observe that a 2 > m 2 /a 2 + 1, see Fig. 10a. In particular, the weaker inequality 
(D6) is satisfied. Thus, for reasonable choices of the parameters e, F n , we find p(0) > 0. 

One may wonder whether p(0) is always positive, irrespective the choice of e, i 7 ^, m n . Recalling that the boundary 
value problem (15) is characterized by m-n only, which follows from substituing x — r e F^ in (15), we observe that the 
inequality (D6) is satisfied at least up to fh^ < 25 (where we start to arrive at the limits of our numerical procedure), 
see Fig. 10b. Our observations suggest that if the boundary value problem (15) has a solution, then it yields p(0) > 
(and, as we explicitly checked, satisfies (10), and yields a negative d-y). It seems that p(0) > is a property of the 
boundary value problem (15), though we are able to verify this only in a finite range of the (positive) parameter m w 
(which, however, goes far beyond what is needed to describe the physical situation in the Skyrme model). 

Another interesting observation is the linear asymptotic behaviour a 2 f» Am, + B with A w 1.688 and B m 0.865, 
which sets in already at moderate values of fh, see Figs. 10a and b. 
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a 2 vs. m 2 /a 2 + l (a) a 2 vs. m 2 % /a 2 + 1 (b) 




0.1 0.2 0.3 0.4 m n 5 10 15 20 m K 



FIG. 10: Solid lines: The dimcnsionlcss quantity a = F'(0)/(eF n ) as function of = m 7r /(ei ; V) in the region of (a) small m-n < 0.5, 
and (b) up to m, n < 25. Dashed lines: The function ml^/ct 2 + 1 which is smaller than a 2 demonstrating that p(0) > 0, cf. Eq. (D6). 



APPENDIX E: CHIRAL PROPERTIES OF THE FORM FACTORS 



In this Appendix we derive the results presented in Eqs. (50-54). 
In particular, we shall see that these properties are determined by 
the long distance behaviour of the profile function (39), i.e. that 
they are of general character. 

For that we notice that it is possible to rewrite the profile func- 
tion as 



F(r) — 2 arctan 



R 2 

/(r, TOjt) (1 + m^r) exp(-m 7r r) 



(El) 




r in fm 



FIG. 11: The function /(r,m w ) defined in Eq. (El). 



with /(r, m^) satisfying the following properties. First, as r — > oo 
the function /(r, m^) — > 1 while all its derivatives (with respect 
to r and/or m T ) tend to zero. This guarantees the correct large 
distance behaviour of F(r) in (39). Second, the soliton equation 
of motion (15) demands that this function has a Taylor expansion 
around r = of the form /(r, m^) = Y] c„ .(m„) r" in which only 
odd powers n = 1, 3, 5, ... appear. I.e., it vanishes for small r as 
f(r, m n ) oc r, which guarantees the correct boundary value for the 
soliton profile -F(O) = it in the Ansatz (El). 

The function f(r, m w ) can be determined uniquely from (El) after Eq. (15) has been solved for F(r) with the 
boundary values ^(0) = tt and F(r) — > as r — > oo. Alternatively, one may insert the Ansatz (El) into Eq. (15) and 
solve directly the new differential equation for f(r, m w ) subject to the new boundary conditions: lim r _>o f( r , Y>hr)/ r = 
const and f(r,m„) — * 1 as r — > oo. We stress that by introducing (El) we merely change the notation but do not 
introduce any approximation. Fig. 11 shows /(r, m^) for m x = and 138 MeV. 

We discuss first the slope of the form factor d\(t) at zero momentum transfer. The general expression for ^(0) in 
terms of s(r) follows from expanding Eq. (5) for small t and reads 



4(0) 



Mn 
42 



d 3 rr 4 s(r) 



M N 
16 



d r r p(r) . 



(E2) 



The second relation in (E2) follows from integrations by parts and the relation (9). From the large distance asymptotics 
of p(r) and s(r) in (42, 43) it is clear that ^(0) diverges in the chiral limit, and is well defined only for ^ (and 
that only for m w ^= manipulations like the integration by parts are well defined). 

We consider first the expression (E2) for d'^O) in terms of s(r). For our purposes the model expression (27) for the 
distribution of shear forces can be conveniently decomposed as s(r) = sa(t) + sb(t) where 



F 2 
4 



F'{rf 



sin 2 F(r) 



s B {r) 



F(r) 



F'{r) 2 



] F{r) 



(E3) 
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In the chiral limit ss{r) oc -jj at large r, and gives a contribution to the slope of di(t) at t — which is well defined 
and finite in the chiral limit. Thus, denoting this contribution by d^(0)s, we note that 

oo 

<(0)b = ~ |d 3 rr 4 SB (r) = 0(m°). (E4) 
o 

In fact, what makes d'i(0) divergent in the chiral limit it is the part sa{t) oc -\ at large r. Let m v ^ in the following. 
We perform the integral over r in a spherical box of radius D, and take D — > oo later. Inserting the expression (E3) 
for s^i(r) with the Ansatz (El) into (E2) and making the legitimate (because m-n ^ 0) substitution r — > a/m n yields 

Dm, 

4^|^ daGl{a ^ } (E5 ) 
o 

where we defined the function (with f'(y,m % ) = -§^f(y,Tn w ) for brevity) 

dia^) = _ x [( a 4 +4a 3 + 7a 2 +6a + 3 ) /( ^ jTO7r) 2 

a 4 + m%R% (1 + a) 2 f(^,m„) 2 e- 2a 

2a (a 3 + 3a 2 + 4a + 2) /(-£-, )/'(;£-, m n ) + a 2 (a + l) 2 /'(^, m,) 2 J- . (E6) 



In Eq. (E5) for 7^ we may safely take D — > 00. This step is trivial in our analytical calculation, however, it is a 
subtle issue for calculations (in models or lattice QCD) performed in a finite volume. The correct chiral behaviour is 
obtained only by first taking the infinite volume limit, and considering then the limit of pion masses becoming small. 
The m^-expansion of the integral (E5) reads 

D OO 

daG 1 (a,m*) = Jda e~ 2a (a 4 + 4a 3 + 7a 2 + 6a + 3) + 0(m n ) = 7 + 0(m w ) , (E7) 


where for the leading term we made use of the fact that f(y,m n ) — > 1 for large arguments y, while its derivatives 
go to zero. Eliminating Ro in favour of qa and according to (39) we obtain the result quoted in (51). We may 
repeat this exercise exploring the relation of d'i(0) in terms of the pressure in Eq. (E2). The calculation is analog to 
that presented above and yields the same result (51). Thus, our procedure to derive the leading non-analytic chiral 
contribution to 0^(0) respects the model equations of motion. 

Next let us discuss the chiral expansion of d\ = d\(t)\t=o- Since s(r) and p(r) oc ^ at large r, see (42, 43), di takes 
a finite value in the chiral limit. The next term in its chiral expansion is linear in m w which can be found as follows. 
We use the relation (11) for d\ in terms of s(r) and the presentation (El) for the soliton profile, and interchange the 
order of differentiating with respect to m w and integrating over r. Substituting r — > a/m^ we obtain 

dd 1 1 f 3 2 ds(r) 4irM N FZR* °[ 

= -- d d rr z — — = J da G 2 (a,TO 7r ) , (E8) 

om^ 6 J om n 6 J 



where G2(a, m T ) is a function similar to (but much lengthier than) that in Eq. (E6), whose small m^-expansion reads 

00 00 

2a / 3 , r>„2 5 



da G 2 (a, m n ) = -2 /da e' 2a (a 3 + 2a 2 + a - l) + 0(m n ) = -- + 0(m n ) , (E9) 

which yields the chiral expansion of d\ quoted in Eq. (50). Starting from the expression (11) which expresses d\ 
in terms of p(r) one reproduces the result in Eq. (50) by a similar calculation. We stress that also the calculations 
leading to (50) respect the model equations of motion. 

Next we consider the mean square radius of the energy density defined in (45). Since Too(r) oc \ at large r, 
see Eq. (41), the mean square radius of the energy density has a well defined chiral limit. Its first chiral correction 
appears already in linear order in m n and the coefficient of this leading non-analytic term can be found by means of 
a calculation which is completely analog to that in Eq. (E8) and yields the result quoted in Eq. (52). 

Finally, we consider the mean square radius of the angular momentum density which is defined in (46). A calculation 
analog to that in Eqs. (E3-E7) yields the result quoted in Eq. (54). 
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APPENDIX F: SCALING RELATIONS 



In our approach — in which we consider 1/N C corrections to some quantity if and only if the respective leading 
order contribution vanishes, see Sec. Ill — the dependence on the model parameters and e is trivial and given by 
the following scaling relations 

*-°(£> *<=-Mt> d '=^ G te> (F1) 

The last result holds for all mean square radii. The functions G(y) are different for the different quantities. 

If one chooses e, differently from the parameter fixing used in this work in Eqs. (20, 21), then the relations (Fl) 
allow to rescale the results for the respective quantities in terms of the new parameters. Strictly speaking one needs 
for that the m^-dependence of the respective quantities which we have provided in Sec. VII for quantities of interest. 
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